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Abstract 



'^ , We consider the exchange statistics of vortices, each of which traps an odd number {N) of 

O I Majorana fermions. We assume that the fermions in a vortex transform in the vector representation 

of the SO{N) group. Exchange of two vortices turns out to be non-Abehan, and the corresponding 
^ . operator is further decomposed into two parts: a part that is essentiahy equivalent to the exchange 

m 

C"^ \ operator of vortices having a single Majorana fermion in each vortex, and a part representing the 

Coxeter group. Similar decomposition was already found in the case with iV = 3, and the result 

m 

^^ ■ shown here is a generalization to the case with an arbitrary odd N . We can obtain the matrix 

representation of the exchange operators in the Hilbert space that is constructed by using Dirac 
-^-. _ fermions non-locally defined by Majorana fermions trapped in separated vortices. We also show 

rS 

H ■ that the decomposition of the exchange operator implies tensor product structure in its matrix 

C^ ■ 

representation. 
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I. INTRODUCTION 



There has been consider ab 
vortices in superconductors 
with non-Abehan statistics 12 



e interest recently in zero-energy fermion modes trapped inside 



ian anyons. They 
[a, |6| and have been 



Vortices in a chiral p-wave superconductor are endowed 
3| because of the zero-energy Majorana fermions inside them 
^. Excitations which obey non-Abehan statistics are caUed non-Abe 
are expected to form the basis of topological quantum computations 
investigated intensively [7|. A recent classification of topological superconductors clarifies 
;he condition that vortices have zero-energy Majorana (or Dirac) fermions in their cores 
8L9|. One remarkable development is the fact that non-Abelian anyons in three dimensions 
10| can be realized by monopoles with Majorana fermions trapped inside their cores, and 



give a new non-Abelian statistics, the projective ribbon permutation statistics 11 |. 

More recently, the non-Abelian statistics of vortices with additional 5*0(3) symmetry 
has been investigated and shown to have a novel structure [12] • In this case, the Majorana 
fermions form the vector representation of the 5*0(3) group. It is shown that the represen- 
tation under the exchange of two vortices is written by the tensor product of two matrices. 
One matrix is identical to the exchange matrix for vortices with a single Majorana fermion in 
each core, found by Ivanov |3| modulo trivial change of basis, and the other matrix is shown 
to be a generator of the Coxeter group, which is a symmetry group of certain polytopes in 
high dimensions [l3|. 

Such vortices can be physically realized in a "color" superconductor. At extremely high 
densities and low temperatures, which could be achieved in the cores of compact stars, 
hadronic matter undergoes phase transition into quark matter, and it is expected to exhibit 
the color superconductivity ij]. Unlike the ordinary superconductivity in a metal, color 
superconductivity is induced by condensation of diquarks, which are pairs of two quarks 
(fermions having "colors" and "flavors"). The order parameter of the color-flavor locked 
(CFL) phase is given by a 3 x 3 matrix, 

^ai = (^ap^,eijk{{q^Y^C-f5 (g)^), (1.1) 

where q is the quark fleld, a, /3, 7 = r, g, b (z, j, k = u,d, s) are color (flavor) indices, and the 
transpose is employed with respect to the spinor index. At asymptotically high densities, 
the ground states are expected to be in the CFL phase, in which the diquark acquires an 
expectation value like $q,j = A6ai, where A is a BCS gap function. The original color 



SU{3)c and flavor SU{3)f symmetry break down to an 5'f/(3)c+F symmetry, the elements 
of which are given by "locked" rotations of color and flavor, $ — t- [/$f/~^. This structure 
is similar to the Balian-Werthamerer (BW) state of the ^He superfluids, in which the order 
parameter is invariant under the locked rotations of spin and orbit states. We can argue 



that there exist topologically (and dynamically) stable vortices [15| in the CFL phase, by 
examining the symmetry-breaking pattern. The vortices in the CFL phase are color flux 
tubes as well as superfluid vortices, since both local and global symmetries are broken 
in the ground state. In the presence of a vortex, the order parameter takes the value 
like $(r) = Adiag{f{r)e^^,g{r),g{r)} where r is the radial coordinate and /(r) and g{r) 
are functions of r. This kind of vortex solution breaks the SU{3)c+f symmetry down 
to S'f/(2)c+F X U{1) symmetry inside the core 16l-l2l| and fermion zero modes belong to 
representations of SU{2)c+f- It has been shown 22|, |23|] that one CFL vortex has triplet and 
singlet Majorana fermions inside it 27(|. Thus vortices in a color superconductor provide 
an example of the system with fermion zero modes in the vector representation of 5*0(3) in 
their cores, since the triplet of SU{2) is equivalent to the vector representation of 5*0(3). 



However, it should be emphasized that the results obtained in Ref. [12| do not depend on 
details of speciflc models. The only assumption adopted there is that a vortex has Majorana 
fermions which transform in the vector representation of 5*0(3), and therefore we expect 
that such a system can be found in condensed matter systems such as exotic superconductors 
or ultra-cold atomic gasses. 



In this article, we generalize the results of Ref. [12[ obtained for 5*0(3) to the case of 
SO{N) where A^ is an arbitrary odd integer A^ > 3. We discuss the exchange statistics 
of two vortices having A^ Majorana fermions trapped in their cores, which are transformed 
in an SO{N) group. Notice that the SO{N) symmetry is the maximum symmetry in the 
presence of A^ Majorana fermions which are real flelds. We discuss the exchange statistics 
in both the operator and representation levels. We flnd that both the operator and the 
matrix representation of the exchange operation generally have factorized structures for 
arbitrary odd A^. In particular, the matrix representation is written as a tensor product 
of two matrices, as previously found in the 50(3) case. One of the two matrices is the 
one discussed in Ref. J3| in the case that a single fermion is trapped in each vortex. We 
show that the other is again a generator of the Coxeter group, as in the case of 5*0(3). 
When one Majorana fermion is topologically protected in the vortex core such as in class- 



D topological superconductors, it is robust against perturbations and remains zero energy 
under the exchange operation. In addition to that, A^ Majorana fermions remain zero energy 
as far as the SO{N) symmetry remains intact. 

This paper is organized as follows. In Sec. [TTl we review the non-Abelian statistics in 
the case of a single Majorana fermion studied in Ref. |3|. In Sec. IIIIl after summarizing 



the previous work on S0{3) case [12|, we present generalization to the case of multiple (A^) 
Majorana fermions with SO{N) symmetry. We explicitly show factorization of the exchange 
operators into the known part similar to the case with a single Majorana fermion, and the 
part corresponding to the Coxeter group. We also show an interesting decomposition of the 
Majorana fermion operators, which clarifies the action of the Coxeter group on Majorana 
fermions. In Sec. |V] we discuss the relation between the decomposition of the exchange 
operator and its matrix representation. Section |Vl] is devoted to summary and discussion. 
In Appendix |A] we prove the decomposition theorem of the exchange operation of vortices. 

II. VORTICES WITH A^ = 1 

We briefly review how non-Abelian statistics emerges for a set of n vortices, each of which 
contains a single Majorana fermion in its core |3|. This provides the simplest example of 
non-Abelian statistics of vortices, but as we will see later we can always identify the same 
structure even for the case with multiple Majorana fermions as far as the number of fermions 
A^ is odd. 



Consider an exchange of two vortices in a system oi n = 2m vortices 28|. Each vortex 
has a single Majorana fermion localized at the core of it, and one can specify the position 
of a vortex on a two-dimensional plane (we label the vortices k = 1, ■ ■ ■ ,n). Notice that 
the trapped Majorana fermion has zero energy, which gives rise to degeneracy of the ground 
(lowest energy) state. Since the existence of Majorana fermions are topologically guaranteed, 
the degeneracy is not disturbed by small perturbations, and hence we treat the exchange of 
vortices as an adiabatic process. 

Let an operation T^ be defined as an exchange of the k-th and {k + l)-th vortices, in 
which the [k + l)-th vortex turns around the k-th vortex in an anti-clockwise way. All the 
exchanges of two vortices are realized by successive application of the exchanges of adjacent 



vortices Tk, and they form a braid group Bn- They indeed satisfy the braid relations: 

TkTiTk = TiTkT, for \k-i\ = l, (2.1) 

TkT, = TiTk for \k-i\>l. (2.2) 

Recall that the vortices are accompanied by Majorana fermions, and one can express the 
action of Tk on Majorana fermions as a transformation. To this end, we define a Majorana 
operator 7^ corresponding to the Majorana fermion in the k-th vortex 3|, satisfying the 
self-conjugate condition (7^)^ = 7^ and the anti-commutation relation {7fc,7£} = 26ki (the 
Clifford algebra). Then, the operation T^ induces the following transformation 



3: 



T.:^^'*"^'*" . (2^3) 

7fc+l -> -7fc 

with the rest '-fe {£ ^ k, k + 1) unchanged. One can explicitly check that the transformation 
(12. 3p satisfies the braid relations (12. ip and (12.20 . The minus sign in the second line is 
essential for the non-Abelian statistics because it gives T^ = 1 (the Bose-Einstein or Fermi- 
Dirac statistics give just T| = 1). The transformation (12.30 is realized by the unitary 
operator 

Tk = exp (^-7fc+i7fc) = 7f (^ + 7fc+i7fe) ■ (2.4) 

Indeed, one finds 

TklkTk^ = 7fc+i , (2.5) 

Tk-fk+lTk^ = -Ik , (2.6) 

rk7irk'=li ii^k,k + l). (2.7) 

We call the operator Tk for this single Majorana case the 'Ivanov operator' since it was first 
found by Ivanov [3|. One can explicitly see that this transformation is indeed non-Abelian 
in the matrix representation of Tk- To construct the Hilbert space on which the operator 



Tk acts, we define Dirac fermions 29[ by using two Majorana fermions at adjacent vortices 
"^K = {'l2K-i + il2K)/'2 with K = 1,- ■ ■ ,171. These Dirac fermion operators satisfy the usual 
anti-commutation relations, 

{^j,, ^[} = 6kl, {"^k, ^l} = {^1,, ^U = 0. (2.8) 



If one defines "^k and \1/)^ as the annihilation and creation operators, respectively, then 
one can construct the Hilbert space by acting the creation operators '^^s on the "Fock 
vacuum-state" |0) defined by ^E'/^lO) = for all K. Within this Hilbert space, the operators 
Tfc's are now expressed as matrices which we call the Ivanov matrices. The Ivanov matrices 
contain off-diagonal elements representing the non-Abelian statistics. Explicit forms of the 
Ivanov matrices for two and four vortices are given in Ref. p|, and also in Ref. [l^] with a 
different expression, both of which are related to each other by unitary transformations. 



III. VORTICES WITH MULTIPLE MAJORANA FERMIONS iV > 3 

Now let us turn to the case with vortices having multiple Majorana fermions in their cores. 
Each vortex traps A^ Majorana fermions 7^ (a = 1, ■ ■ ■ , A^) which transform in the vector 
representation of SO{N) symmetry, hence we call them non-Abelian Majorana fermions. 
We consider A^ to be an arbitrary odd number, and this is a generalization of the simplest 



non-trivial case A^ = 3 in Ref. 



12|. 



The non-Abelian Majorana operators 7^ satisfy the self-conjugate conditions and the 
anti-commutation relations: 

iiky=ii {iiii} = '^5^%i- (3.1) 

We define the exchange of the k-th and {k + l)-th vortices so that the Majorana fermion 
operator with each a transforms in the same way as the case with a single Majorana fermion 



(see Eq. (El) |3Q|: 



Tfc : <( ^' ^ ^'+' for all a , (3.2) 

7fc+i ^ -it 

with the rest 7° {£ ^ k, k + 1) unchanged. The exchange operations T^'s satisfy the braid 
relations (12. ip and (12. 2p . Their action on 7^'s can be represented in terms of non-Abelian 
Majorana operators 7^'s in the following way. Since the transformation (13. 2 p for each a 
is equivalent to the single Majorana case, one can use the same expression for the unitary 
operator which induces the transformation for each a: 

Tl ^ exp i^-iU.ii) = ^ (1 + 7^i7.1 • (3.3) 
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Thus, the exchange operator for the vortices having multiple fermions should be represented 
as the product of them: 

N 

r^-^-^ (3-4) 

a=l 

This exchange operator is SO{N) invariant as shown in the next section. One can check 
that the operator r], applied to 7" indeed generates the desired transformation fl3.2p : 

rrht{rrY' = lU., (3.5) 

rf 7.V(rf ')- = -7.% (3-6) 

rrhi{rrr'=li {i^k,k + l). (3.7) 

This transformation is again non-Abelian, which is explicitly seen in the matrix representa- 
tion. 

To obtain the matrix representation, we can perform the same procedure as in the case 
with single Majorana fermions. Namely, by defining the Dirac fermion operators 

n = ^(72^-1 + ^12k), < = ^(72^-1 - ^12k) (3.8) 

which satisfy {K,L = 1,- ■ ■ ,m) 

{vl/^,vl/^t} = 5^z.5'^^ {vl/^,vl/^} = {M/^t ,]>6t|^Q^ (3 9) 

we can construct the Hilbert space. Then, we can find matrix representation of r^ . In 



Ref. 



12| . three of us explicitly constructed the Hilbert space for the case of A^ = 3 and 



[si 
n = 2,4, and found matrix expression of rj, according to different representations of 5*0(3). 

The matrices have off-diagonal elements and thus they are non-Abelian. In principle, one 

can do the same thing for an arbitrary odd A^. In the present paper, however, we first 

look into interesting structure of the operator r^, , which was also suggested in the previous 



JJV] 



paper 121]. Namely, the operator r]. itself can be decomposed into two parts. Then, we 



discuss the relation between the decomposition of r^ and the matrix representation of r^ 



m 



Sec.lYl 



IV. THE COXETER GROUP FOR MULTIPLE MAJORANA FERMIONS iV > 3 



In the previous analysis [12| with A^ = 3, we found that the matrix representation of rj^ 



[3] 



in the four- vortex sector can be decomposed into a tensor product of two matrices: one is the 



same as the Ivanov matrix for the single Majorana fermion case, and the other is identified 
with generators of the Coxeter group. We also found that the similar decomposition is 



possible at the operator level [12|. Namely, one can express r^ as a product of two distinct 
operators which give rise to the corresponding matrices in the matrix representation. In 
this section, we discuss that such a decomposition at the operator level holds even for an 
arbitrary odd number of A^. 

Before we go into details, let us define some useful notations. We first define a composite 
operator F^ by 

n^ii+iii. (4.1) 

which have the following properties, 

{VU\} = -25ki for a = 6 and |A;-/| < 1, (4.2) 

[r^,rf] =0 for a ^6, or (a = 6 and |A; - /| > 1) . (4.3) 

For later convenience, we define another composite operator (1 < n < A^) 



pV"; — 

'- k — 



/^n2 ^ ^ ^ "b„+i ''bjv 7fc+l Ik+llk Ik ' \^-^) 



where e"^'"°^ is the completely antisymmetric tensor. It is evident that the operators F^*^ s 
are SO{N) invariants for all n. For example, for A^ = 3, 

F« = F1+F^ + f3, Ff = FiJ^ + r^F^ + F^F^, T^^ = TlTlTl (4.5) 

With those composite operators, the exchange operator rj. can be represented as 



^k 



[TV] 



^)n(-r..e.p{^Erf (;^)grr. («) 



In the final expression, we confirm that the operator r], is S'0(A^)-invariant. 

A. The case of 50(3) 

The simplest nontrivial case with A^ = 3 provides us with useful information which is 
helpful in discussing the decomposition for the general case A^. Let us first recall that 
the operator r]. ~ has been found to be decomposed into two parts (see Appendix B in 



Ref. 



12|): 

rP = a^'hf , (4.7) 



where both of the operators (t[, and hf. are given in terms of the Majorana operators 7^ as 



[3] 



- (1 - il+iil+iihl - il+iil+iihl - il+iil+iihi) , (4-8) 



and 



4'] = ^ (1 - lUiiUiil+iihhi) ■ (4.9) 

One can also rewrite them compactly in new notations as (see Eqs. fl4.4p and f l4.5p ) 

4" = i(i + rf)- "'' = 7i (1 + r*") ■ '"•" 

[31 [31 [31 

Note that a^ and h\ are commutative for any pair of k and L Thus, r^ is the product of 
a\^ and h)^ . It was shown that the operators cr[ s satisfy the Coxeter relations 3l|, 

{^ff = 1, (4.11) 

[afaff = l for \k - E\ = I, (4.12) 

{afaff = 1 for \k - i\ > 1. (4.13) 

The Coxeter group found here for the n = 2m vortices is classified as A2m-i- It is also known 

as the symmetric group S2m, which is the symmetry group of a regular (2,771 — l)-simplex. 

[31 
In contrast, the other part h^, works in the same way as the Ivanov operator r^ defined 

in Eq. (12. 4p . although /ij, has a more complicated structure. This is naturally understood 

if one notices that the operator jI^I'jI is 5*0(3) invariant (singlet), and thus it treats three 

Majorana fermions together as if it does not have any 5*0(3) index. This picture is very 

useful when we consider the general case with A^. 

B. The case of SO{N) with arbitrary odd N 

The previous analysis suggests that it is possible to decompose the full exchange oper- 
ator tI into two parts. This is indeed the case. We find that the operator r^ can be 
decomposed as 

rr'^arAfl, (4.14) 



where cr[ and h\. are defined by using the notation introduced before as 



^i + n^^ + rr + --- + rr^M, (4.15) 






Note that a[, and h\, are SO{N) invariant, and a\. and /i^ are commutative for any k 
and L One can readily verify the decomposition f l4.14p . By using Eqs. f l4.15p and fl4.16p . 
one can check that the product a\. h\. is equal to the last equation in Eq. (14. 6 P if one uses 
Eqs. (??) and (lilall . 

First of all, let us discuss the properties of /i[. . The analysis for the case A^ = 3 suggests 
that if one treats multiple Majorana fermions in a vortex in a unit, then the action of t\, ' 
will be essentially equivalent to the Ivanov operator. This motivates us to introduce the 
following "singlet Majorana operator" locally defined on the /c-th vortex: 

Ik = ^e^^(^-i)6-"™7r7r ■ --iT , (4.17) 

which is manifestly invariant under the SO{N) transformation. The phase factor is included 
so that the operator becomes self-conjugate (7^)^ = 7^, and satisfies the Clifford algebra 
{lk:7i} = 25fc^. Notice that these properties of 7^ are the same as those of a single Majorana 
operator. For A^ = 3, one finds 7^ = i'jl'jl'jl, and the operator hj^ can be compactly 
expressed as h[. = y={l +7^+17^), which has the same structure as the Ivanov operator 
(12. 4p . For arbitrary odd A^, we find that h[. can also be expressed as 

'71 \ 1 



/if] = exp (^7.+i7.) = 7f (1 + 7.+17.) , (4.18) 



V2 
by noting the relation F^ = 7^+17^. Then, h[. works on 7^ as 

C^7.(/^f')-^=7.+i, (4.19) 

4'"7.+i(4''y^ = -Ik, (4.20) 

/^r7.(Cy^=7. {i^k,k + l), (4.21) 

for an arbitrary odd A^ > 3. Interestingly, this is the same as the transformation (12. 5p - 
(12. 7p induced by the Ivanov operator with A^ = 1. Therefore, the operator h\. for the 
singlet Majorana operator 7^, is equivalent to the Ivanov operator r^ for the single Majorana 
operator 7^.. 

10 



Next, we discuss the properties of the operators a\. . Similarly to the A^ = 3 case, a\^ are 
generators of the Coxeter group. Indeed, it can be confirmed that cr[ satisfy (see Appendix 
\M for details) 

(4.22) 
(4.23) 
(4.24) 



{^f'f - 


= 1, 








af af )3 = 


= 1 


for 


k- 


-i =1 


4%r^r - 


= 1 


for 


k- 


-i > 1 



and these relations induce the same Coxeter group for n = 2m vortices. We thus have found 
that, for an arbitrary odd A^, the operators a\. s again obey the Coxeter relations of A2m-i- 

In fact, one can "derive" the decomposition (I4.14p by first defining h\. in terms of the 
singlet Majorana operators 7;, as in Eq. (I4.18p . and then assuming the factorized form (I4.14p . 
The operator a^ = tI {h\. )~^ thus obtained indeed coincides with Eq. (I4.15p . and the 
commutativity of a^ and h[ for any k and i confirms the assumption of factorization. 

We thus have found that, for an arbitrary odd A^, the exchange operator r^ is expressed 
as a product of a generator of the Coxeter group A2m-i (for the vortex number n = 2m) 
and the Ivanov operator for a single Majorana fermion. 

C. Decomposition of the Majorana operators 

Let us recall that the exchange of the Majorana operators 7"'s is originally defined as 
the operation T^ in Eq. (13. 2p . It is not apparently clear how the decomposed structure 
of the operator t\. indeed works in the exchange operation T^. To understand this, it is 
instructive to notice that the Majorana operator 7^ can be rewritten as 

ll = ltlk. (4.25) 

where 7^ is a composite operator in the vector representation of SO{N) defined locally on 
the k-th vortex as 

Ik = ^^^^e^^(^-i)6----7r ■ ■ ■7r-^ (4-26) 

and % is the singlet Majorana operator defined in Eq. (I4.17p . The two operators 7^ and 7^ 
commute with each other for any pair of k and i. This expression allows us to extract, from 
the Majorana operator 7^, the part of a singlet Majorana fermion 7^ whose properties are 

11 



well understood. Notice that 7^. (7^) is composed of an odd number A^ (an even number 
A^ — 1) of Majorana fermion operators. 

Since 7^ and 7^ are expressed in terms of the original Majorana operator 7^, one can 
immediately find how they are transformed in the exchange T^. Namely, we find the trans- 
formation of 7^ and 7;, by T^ as 



7fc+i ^ 7fc 



Tk:{l^ ^ 2''^\ for all a, (4.27) 

without a minus sign, and 

Tfc : <( ' " -^^ ^ fo^ ^YY ^ ^ (428) 




with a minus sign, while the rest 7^ and 7^ (£ 7^ A; and k + 1) are unchanged 32|. It is easily 
checked that the simultaneous transformation of 7^ and 7;. reproduces the transformation 
of 7^ in Eq. (13. 2p . Therefore, we observe from Eq. (I4.27P that 7^'s are transformed by a 
symmetric group S2mi or the Coxeter group of A2m-i (for the vortex number n = 2m), and 
from Eq. (I4.28P that 7^'s are transformed as in the same way as the Ivanov operators for 



single Majorana fermions 331 ]. 



The exchange properties of 7^ and 7^ in Eqs. (I4.27P and (I4.28P can be discussed at the 
operator level. Because a^ and 7; {h[, and 7^ ) are commutative for any pair of k and i, 

[4'",7.] = [/^r,7.1 = 0, (4.29) 

the transformation r], 7"(t], )~^ is decomposed as 



rri^rrr = |^r^7."«r^| |/^r^7.(/^rr^) ■ (4.30) 

Therefore, 7f and 7; are transformed by a], and h[, , respectively. From Eqs. (I4.15P and 
(I4.26p . 7^ is transformed as 

-r'iti-rr=iu. (4.31) 

-r7,\,(-r)-^ = 7^ (4.32) 

af7r(4'")-' = 7," {i^k,k + l), (4.33) 

without a minus sign. Thus, the operator cr[. acting on 7^" reproduces the transformation 
(I4.27p . We note that a^ can be expressed in terms of 7^ only. On the other hand, 7;, 

12 



is transformed by the operator h[. like a single Majorana fermion as demonstrated in 
Eqs. f l4.19p - p.2ip . and hence h[. reproduces the transformation f l4.28p . 

To summarize this subsection, in correspondence to the product of r], = cr[ /ij, , the 
Majorana operator 7^ is also expressed by the product of the two parts, 7^ obeying the 
Coxeter group given by a^ and ^f. obeying Ivanov's exchange given hy h[ . 



V. OPERATOR DECOMPOSITION AND MATRIX REPRESENTATION 

So far, we have been discussing the factorized structure of the exchange operation of two 
vortices at the operator level. Everything was written in terms of the Majorana operators, 
and the decomposition into the Coxeter and Ivanov parts was naturally understood by using 
the Majorana operators. In contrast, in order to obtain the matrix representation, the usual 
procedure is to define Dirac fermion operators and use them in constructing the Hilbert 
space. Since the Dirac fermion operators are constructed from two Majorana fermions 
located separately at different vortices, it is not trivial if the factorized structure at the 
operator level is preserved in the matrix representation. For example, the Dirac fermion 
operator defined in Eq. (13. 8 p can not be decomposed similarly as the Majorana fermion 
operator as shown in Eq. f l4.25p . In this section, we are going to show that the decomposition 
holds even in the matrix representation in a suitable basis, and discuss the relationship 
between the the decompositions in the operator- and matrix-representation levels. 

In the following, we discuss the case with A^ = 3 for simplicity, but present the procedures 
to obtain the matrix representation so that they can be easily extended to the case with any 
oddiV. 

A. Construction of the Hilbert space 

Let us consider an even number n = 2m of vortices. Then we can construct m Dirac 
fermion operators \1/^ (a = 1,2,3; K = 1,- ■ ■ ,m), given in Eq. (13. Sp . in the vector repre- 
sentation of 5*0(3). The Fock vacuum |0) is defined by the Dirac fermion operators as 

^^|0) = , for all K anda. (5.1) 



One can construct the basis of the Hilbert space by acting the Dirac fermion 



operators \1/;^ 



at 



on the Fock vacuum |0). The explicit forms of the basis were given in Ref. 12|] for A^ = 3 
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and n = 2,4. Now let us construct the Hilbert space in a way different from Ref. [12|. To 
this end, we define the number operator for the triplet Dirac fermions of the i^-th pair of 
vortices by 

M^ = ^i^K- (5-2) 

A generic state for the K-th pair of vortices can be expressed in terms of the eigenvalues of 
this number operator as 

Wh^hJ^K)K (5.3) 

where M^ = or 1 is the occupation number of the fermion created by \1/^ (here we use the 
same character for the operator and its eigenvalues). Then the basis of the whole Hilbert 
space is composed of the tensor product of the states for each K, 

|(8)|A/'i,Arl,Arl)^|. (5.4) 



K=l 



B. Singlet Dirac operators 

When A^ = 3, the singlet Majorana operators 7^. defined in Eq. (I4.17p are given by 

Ik = ^e'^'^lHlt ■ (5.5) 

By using these operators, we define singlet Dirac operators 

'^K = - {12K-1 + mx) , ^K=2 (72X-1 - mx) ■ (5.6) 

These operators play a particular role when they act on the states. Let us examine the 
action of the singlet Dirac operators "^k^s and ^I/j^'s on the Hilbert space (15. 4p . which is 
constructed by acting the triplet Dirac operators \l/^'s and \l/^'s on the Fock vacuum |0). 
To this end, we express the singlet Dirac operators "^k and "^j^ in terms of the triplet Dirac 
operators \I^^' and \&^. For that purpose, we note that the singlet Majorana operators ^2K-i 
and ^2K defined in Eq. (I4.17P can be written in two ways: 

72K-1 =^K+% 






14 



(5.7) 



respectively. From these two relations, we can express the singlet Dirac operators ^^ and 
\1/^ in terms of the triplet Dirac operators \1/^ and \1/^ as 



% = -^^e"'' {^K^K^K + ^"ll^K^i + ^X^X^^ + ^K^K^k) ■ (5.8) 



For the later purpose, we define the fermion number operator for the K-th pair of vortices 
in terms of the singlet Dirac operators \1/^ and "^k'- 

Wk = %'^K- (5.9) 

The meaning of this number operator will be clarified below. 

Now we can discuss how the singlet Dirac operators "^k and \1/^ act on the Fock states 
(15. 3 p defined by the triplet Dirac operators \1/^ and \1/^. The action of the singlet Dirac 
operators "^k and \&;^ on \Af]^,M^,Af^)K can be read off from Eqs. (15. 8 P as 

^^|0,0,0)i, = |l,l,l)i^, 

^i^|l, 0, 0)^ = ^,,10, 1, 0)/^ = ^x|0, 0, 1)/^ = 0, 

^^|0, 1, 1)k = |1, 0, 0),^, ^^|1, 0, 1)k = |0, 1, 0),,, ^x|l, 1, 0),^ = |0, 0, 1)k, 

^^|1,1,1),< = 0, (5.10) 

for "ifx, and 

^;,|o,o,o)a' = o, 

^1,10, 0, 1)k = |1, 1, 0)k, %\0, 1, 0)^ = |1, 0, 1),^, %\1, 0, 0),, = |0, 1, 1)^, 

^;,io, 1, 1),^ = ^;,ii, 0, 1),^ = ^;,ii, i, o),^ = o, 

^;^|1,1,1)k = |0,0,0),^, (5.11) 

for "^1^. 

Let us define the total fermion number operator for the J'i'-th pair of vortices as the sum 
of the number operator Af^ in Eq. (15. 2p of the triplet Dirac fermions: 

N 

J'k^Y.^k- (5.12) 

a=l 
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Then we can deduce that the states annihilated by "^k are those with odd eigenvalues of 
J^K (see the second and fourth lines of Eq. fIS.lOp ). while those with even eigenvalues of Tk 
are annihilated by \E'^ (see the first and third lines of Eq. f IS.lip ): 

*'J^|even)i^ = 0, '^K\odd)K = 0, (5.13) 

where |even)x and |odd)i^ are eigenstates of the parity operator 

Vk = i-lf", (5.14) 

namely, 

Pif |even)/^ = +l|even)i^, VK\odd)K = — l|odd)i<'. (5.15) 

It should also be noted that, when \E'^ and "^k do not annihilate the state, they create a 
state with opposite parity: 

^j^\odd)K = |even)i^, \E'x|even)x = |odd)j<'. (5.16) 

These facts imply that the parity operator Vk anti-commutes with "^k and \E'^: 

{Vk,^k} = 0, {Vk,%} = 0. (5.17) 

We claim that the fermion number operator Wk is related with the parity operator Vk 

as 

VK = l-JfK. (5.18) 

Namely, the fermion number operator 77 k expresses the parity of the total fermion number 
for each index K. This relation results from Eq. ( I5.13p . which is obvious from the structure 
of the operators shown in Eqs. (15. 8p . 

By repeating the same argument, the above relation can be generalized to an arbitrary 
odd A^ as 

p^ = l_Ar^ for N = M + 3, (5.19) 

Vk=J7k for N = U + l. (5.20) 

C. The tensor product structure of the Hilbert space 

From Eqs. fl5.10p and (15. lip , we can say that the action of the singlet Dirac op- 
erators ^ K or ^ j^ is a kind of "NOT" operation, which "flips" the fermion number 
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for each component of the 5*0(3) vector. Namely, under the action of "^k, one finds 
Wk^-^k^-^k) A' — ^ |1 — -^K' 1 ~-^K' 1 ~ ■^k)k- W^ ^^^ divide the states into pairs, in 
each of which the two states are related by the NOT operation. When A^ = 3, there are four 
pairs: 

|0,0,0)i^o|l,l,l)i^, 10,0,1),^^ |l,l,0)i^, 

|0,l,l)i,o 11,0,0);,, |l,0,l)i,^|0,l,0)i,. (5.21) 

It is important to note that the singlet Dirac operators '^ k and ^ j^ induce the transi- 
tion only within these pairs when they act on states (e.g. a transition between |0,0,0)i, 
and |1, 1, 1)^-), but they never induce an inter-pair transition. This fact is essential in the 
following discussion. 

Now we are ready to discuss how the decomposition of the exchange operator r]. results 
in the tensor-product structure in a matrix representation. From the analysis above, we 
can take the basis of the Hilbert space which is labeled by the parity Vk of the number 
of fermions with index K, and an additional index ttlk which labels the choice of a pair in 
Eq. fl5.2ip . We denote the states by 

[Pk, mK)K = \'Pk)k ® |"^x)x- (5.22) 

Let us consider the matrix elements k^Ki ^^kvu I'P'k^ '^'k)k of the exchange operator r|. ' 
in this basis. Now we show that these matrix elements can be written as the tensor product 
of the Ivanov matrix and the Coxeter matrix: 

K{VK,rnK\Tf \V'K,mj^)K = KirriKWk \m'j^)K ■ K{VK\hk \P'k)k- (5.23) 

Namely, the Coxeter operator a\. acts only on |?7i/<)/,, while the Ivanov operator h\. acts 
only on \Vk)k- To prove this, we show the following two statements for any k and K: 

(i) h\, acts as an identity on \niK)K- 

(ii) a\, acts as an identity on \Vk)k- 

First, we can see that the statement (i) is true in the following way. The Ivanov operator 
h\. is written by the singlet Majorana operators 7^ (see Eq. 04.181) ). and hence by the singlet 
Dirac operators ^ l and \E'^, see Eq. (15.61) . When the singlet Dirac operators \E'l or \&^ act 
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on states, the only state that can be created is the one in which the fermion numbers are 
flipped. So, the action of the singlet Dirac operators ^l or "^^^ never induce an inter-pair 
transition. Therefore, the Ivanov operator h[, does not change the index rriK, which labels 
the pair in Eq. fl5.2ip . 

Next, the statement (ii) is equivalent to the statement that the Coxeter operator a^ 
and the total fermion number operator J^fx commutes. The total fermion number operator 
T^fx is written by the singlet Majorana operators ^2K-i ^^^ l2K^ which commute with the 
Coxeter operator a\. as in Eq. fl4.29p . We thus find that the Coxeter operator a\, and the 
total fermion number operator JJk also commute. 

From these discussions, we have shown that the factorization of the exchange operator 
r|. of vortices into the Ivanov operator h\. and the Coxeter operator a\. results in the 
tensor-product structure in the matrix representation. 



VI. SUMMARY AND DISCUSSION 

We have considered non-Abelian statistics of vortices, each of which has A^ Majorana 
zero-energy states inside its core on which an SO{N) symmetry acts. We have investigated 
how the degenerate states induced by zero modes are transformed under an exchange of 
neighboring vortices. We have shown that, for an arbitrary odd N, the exchange operator 
tI defined in Eq. (13.41) . generating the exchange of two neighboring vortices, can be fac- 
torized into two parts r]. = cr}. /ij, as seen in Eq. fl4.14p . The part which is given by 
h[. defined in Eq. (14.160 is essentially equivalent to the exchange operator introduced by 
Ivanov. If it is expressed in terms of the composite singlet Majorana operator 7;, defined in 
Eq. ( ]4.17p . then it has the same form as the exchange operator r^ in the case of the single 
Majorana fermion. The other operator a^ defined in Eq. (I4.15P is a new part. We have 
proven in Appendix [A] that they constitute the Coxeter group of the type A2m-i (the sym- 
metric group 5*2^) for n = 2m vortices. We have also shown in Sec. |V]that the factorization 
of the exchange operators results in the tensor-product structure in its matrix representation 
in a suitable basis. 

The SO{N) symmetry considered in this paper is the largest symmetry group in the 
presence of N Majorana fermions. Whether the symmetry is SO{N) or its subgroups de- 
pends on the details of the systems. For instance, a higher (pseudo-)spin S representation of 
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so (3) contains 25*+ 1 Majorana fermions, but the symmetry acting on them does not have 
to be 5*0(25' + 1). Also the representation does not have to be irreducible; for instance four 
Majorana fermions may be decomposed into one singlet and one triplet of 5*0(3), but the 
symmetry group does not have to be 50(4). It will be interesting to extend our results to 
general representations including reducible representations. An extension to general groups 
also remains as an interesting future problem to be explored. 

For vortices with an even number N of Majorana fermions with the SO{N) symmetry, 
we have not found any meaningful factorization of the exchange operator r], so far. It 
remains as a future problem to identify the non-Abelian statistics for even A^ Majorana 
fermions. When the symmetry group inside the vortex core is restricted to the unitary 
subgroup U{N/2) C SO{N), N Majorana fermions can be rearranged into A^/2 complex 
Dirac fermions in each vortex. In this case, the situation is rather different because Dirac 
fermions are locally defined, and we do not need to define Dirac fermions non-locally by 
using two spatially separated vortices. Nevertheless we found a rather different kind of non- 
Abelian statistics in the case of A^ = 2 with the U{1) symmetry 25| and A^ = 4 with the 
U{2) symmetry |26|. Identifying the statistics for the general even A^ of A^/2 Dirac fermions 
also remains as a future problem. 

Finally, it will be important to look for actual condensed matter systems realizing multiple 
Majorana fermions in the vortex core, and to study an impact of our results on applications 
to topological quantum computations. 
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Appendix A: Proof of Coxeter relations for arbitrary odd A^ 

In this appendix we give a proof that the operators a\, defined in Eq. fl4.15p satisfy the 
Coxeter relations in Eqs. fl4.23p and fl4.24p . 

1. (aflaf)3 = lfor|fc-£| = l 

In this subsection we show Eq. f l4.23p . Let us note that the cube of the product of r], ' 
and t\,_^^ is written as 

(rf 'rffi)' = (4"'4"VKi ft^i )'. (Ai) 

Since a\, commutes with h)^ and h\_^^, the relation above is written as 

(rrVS.f = (4'"4:\)'(''r'45,)'- (A2) 

The left hand side of Eq. (IA2p is equal to —1, which can be shown as follows. 

N . .3 

^ 1 1 



(A3) 



a=l 



(A4) 



= (-1. 
= -1, 
where in the third line we have used the relation, 

1^ (1 + r^ + Tl^, + TlTl_,,) I = 1 (-1 + r^ + Tl^, + TlTt^,) , 

which follows from the anticommuting property of F^ and r^_,_^. In the fourth line, we have 
again used the anticommuting property. 

On the other hand, we can also show that {hj^ hk_^i)^ is equal to —1, as 

crC)' = {i (i + rr + C, + rf rr,i) }' 

= i (-1 + rr + rS H- rrrw ) i (i + rf + rffi + rrrKl) (^5) 



-1. 
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where we have used 



{i (i + rr + r^ + rf C.) }' = i (-1 + r™ + Tit\ + rrrffi) . (A6) 

From Eqs. (jMl), (M, and (jMD, we can conclude (4^VJ^^)^ = 1 for |A; - £| = 1. 



2. (af]af])2 = lfor |A;-£|>1 

To prove Eq. (I4.24p . we first note that by squaring both sides of Eq. ( I4.14p and using the 
relation 

1 ,. „.A' 1 



— (l + r^)J =-{l + 2n-l) = Tl (A7) 



one finds 

.(TV) _ / [TV]x2 riN) 



n" = K'T rr^ (A8) 



Then, by multiplying (F^ ) ^ on both sides from right, we obtain 

JNU2 



Kr = 1- (A9) 

It follows that (a^. a\ Y = 1 ^^ \k — ^\> 1 since a\. and a\ with |/c — £| > 1 commute. 
We thus have shown that the operators a\. satisfy the Coxeter relation of the type A2m-i- 
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The singlet Major ana fermion found in 
consequently to be non-normalizable 23|] 



22| was in fact shown to diverge at the origin and 



The structure of the Hilbert space in the presence of an odd number vortices is recently 
discussed in Ref. [2J]. 



We use k,i = 1,- ■ ■ , n to label the vortices and the trapped Majorana fermions, while K,L = 
1, • • • , n/2 = 771 to label the Dirac fermions which are constructed from two Majorana fermions. 
One may allow for mixture of indices a under the exchange of two vortices, but here we discuss 
the simplest case where such mixing does not occur. 
A Coxeter group S is defined as a group with distinct generators Sj G S (i = 1,2,3, •••) 



satisfying the following two conditions [l^: (a) s'f = 1 and (b) (s, Sj)"^*-^ = 1 with a positive 
integer rriij > 2 for i ^ j. It should be noted that condition (a) gives rrii^i = 1. Elements rriij 
can be summarized as the Coxeter matrix {M)ij = TUij. In our case of the n = 2?7i vortices, 
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the Coxeter matrix is given by a {2m — 1) x (2m — 1) matrix whose elements are 1 (diagonal 
elements, rrii^i = 1), 3 (adjacent elements, rrii^i^i = rui^i^i = 3) and 2 (all the others). 

[32] In general, a composite operator made by an even (odd) number of the Majorana operators 
is transformed as in Eq. ([071) (Eq. ([08|) ). 

[33] If one considers the case when N is an even number, one may define a composite operator 
by ^'f, = 77Te**^7^.7| • • '7^^- The operator is self-conjugate (7';.)^ = %, but does not satisfy 
the Clifford algebra. Furthermore, since N is even, the operation T^ gives the transformation, 
7'^ — )• ^f._^i, 7fc+i — ^ 7fc and the rest ^'^ {£ ^ k and k + 1) unchanged. Hence the composite 
operator Yfc for even N does not transform like a singlet Majorana fermion operator. 
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